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1. INTRODUCTION 

The functioning of biomolecules has been a field of very extensive research in modern 
biology. One of the interesting questions e.g. is the understanding of transcription and 
copying of DNA molecules. In this contribution we apply recent results we obtained for 
vibrational energy localization in (arbitrary) nonlinear Hamiltonian lattices to models 
which might reflect parts of the dynamics of H-bonds in DNA molecules. We show using 
an abstract level of analysis that there exist generic features of whole classes of models 
with respect to the existence of (very) localized vibrational excitations. By analyzing the 
properties of energy (heat) flows in the presence of these localized vibrations we show 
the localized vibrations to control the heat flow, and thus can define functioning on the 
level of these physical models. Let us however at the beginning spend a few paragraphs 
in order to express our views on the advantages and limitations of theoretical physics 
approaches to the initial problem of functioning in biomolecules. 

The formulation of the problem of functioning of biomolecules (let us for clarity con- 
strain ourselves to DNA transcription) in terms of theoretical physics amounts to describ- 
ing specific nonequilibrium processes of a complex system of many constituents ('parti- 
cles'). These processes might correspond to certain (nonlinear) dynamic processes within 
a choosen physical model. Still measurements of the biological system might reveal cer- 
tain equilibrium properties as well (whatever their connection to the functioning is). A 
good example is denaturation of DNA, which reminds us of a first order phase transition. 
Thus we see two approaches in which abstract theory of many-particle systems might help 
to understand DNA functioning. Both of them start with certain model classes known 
to be good abstractions of different reaUzations of condensed matter (here the crucial 
point is the choice of 'reductionism' of a certain biological constituent to a 'particle', e.g. 
assuming that a base of a DNA-strand is a rigid body, connected to its neighbours by 
certain bonds). The first approach leads us to the definition of necessary conditions our 
models have to have in order to say reveal certain equilibrium properties of DNA. A good 
example is the work of M. Peyrard and coworkers [1,2], who have studied the denaturation 



of DNA. After the answer is found (with respect to the necessary conditions) we will have 
some restrictions on the choice of models and parameters. 

The second approach starts again with the choosen model classes (or it might also start 
from the subclass of models obtained from the analysis using the first approach) and poses 
the question: 'what are the typical nonhnear excitations my class of models posess, and 
how can the presence of these excitations in the original DNA molecule contribute to 
functioning?'. Let us emphasize that the two approaches are not alternative - in fact after 
obtaining the results from the first approach one has to proceed with the second approach 
anyway. In other words, the first approach yields certain constraints (with respect to the 
choice of a model) which have to be obeyed in carrying out the second approach. In 
order to complete the second approach (assuming that we found some typical or generic 
nonlinear excitations or processes) a nontrivial feedback with experimental knowledge of 
functioning of DNA is needed. 

In fact there exists a third way, which is solely computer simulation of models closer 
and closer to the original DNA. Besides the always present computer-assisted limitations, 
this approach is fundamentally different from the first two described ways, because it 
does not pretend to understand a physical phenomenon and then to find its trace in 
the biological system (this will be simply impossible because of the complexity of the 
biomolecule and its functioning). Still this approach might be as (or even more) successful 
in answering questions about functioning of say DNA than the first two ways. 

In the following we will present our results following the second approach, and of 
course we will have to stop at the end because of lack of the mentioned feedback. Still we 
hope that our results and thoughts will help to interprete experimental data and to focus 
and direct future experimental investigations. 



2. H-BOND DYNAMICS 

2.1 Model classes 

As indicated above, we will focus on the dynamics of H-bonds between two bases of 
a basepair in DNA. The simplest approximation of this problem might be to consider the 
bases as rigid masses, connected with their second partner of a basepair through an inter- 
action potential V{z) which provides bound states, and connected to their neighbour bases 
along the strand through a second interaction potential ^{z). We assume discrete trans- 
lational invariance of the system (which will not be crucial, but demonstrates the novelty 
of the results in the best way), i.e. we exclude from our considerations differentiation 
between different bases and/or disorder. For the moment we consider nearest neighbour 
interaction along the strand and assign to every rigid mass one degree of freedom (5±i,i 
where ±1 indentifies the two different strands and I counts the bases along the strand 
{Q could be a rotation angle or just a scalar displacement). The interaction potentials 
are assumed to be given by differences of the involved degrees of freedom {Q+i^i — Q-i,i) 
and {Q±i,i — Q±i,i±i)- Treating the rigid masses classically we use a transformation of 
the original base coordinates and separate the Hamiltonian H into a part Hg which de- 
scribes solely the dynamics of the sum coordinate Yi = Q+i i + Q-i i and a part Hd which 



describes the dynamics of the difference coordinate Xi = Q+i,; — Q-i,; (for details see [1]): 
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The two interaction potentials are given through Taylor expansions around the stable 
groundstate of the system: 

Viz) = y: h'^nz'' , ^{z) = y: 3^"^" • (3) 

n=2,3,... "■• n=2,3,... "- 

We have to choose a potential V{z) which should describe the H-bond energy. The 
physical requirement is that V{z) has a bound state (minimum in z) and unbounded 
states (both because we can break an H-bond by using a finite amount of energy and 
because of the fact of the denaturation transition). Consequently V{z) will have an 
energy barrier which separates the bound states from the unbound ones. If the function 
V{z) is infinitely often differentiable, then it follows that the frequency of oscillation of 
an imaginary particle in this potential will equal y^ for infinitely small energies and 
decrease to zero (continuously) as the energy is increased up to the barrier height. 

If we would consider only small amplitude oscillations of our variables Xi around 
the groundstate X; = of the system, a linearization of the equations of motion Xi = 
—dHd/dXi yields phonon lattice waves as exact solutions (within the hnearization) with 
the dispersion relation u^ = V2 + 4^2sin^(Q'/2), < q < tt. Here q is the phonon wave 
number. Thus we find an optical phonon band with the property V2 < uj"^ < (^2 + 4:4>2)- 



2.2 Nonlinear localized excitations 

These ingredients are sufficient for us in order to make the following statements. System 
(2) supports sets of one-parameter famihes of time-periodic nonlinear localized excitations 
(NLEs) of the form 

Xi{t) = Xi{t + 27r/ui) , Xi^±^ -^0 , 0<ui<y^ . (4) 

The parameter of each family is the fundamental frequency Ui. Time-periodic NLEs 
(4) are exact solutions of the equations of motion if the following condition is met: the 
nonresonance condition kui 7^ Uq has to be fulfilled (here k is an arbitrary integer and 
ujq is an arbitrary phonon frequency). This condition can be certainly fulfilled if (f>2 <^ 
3/4i'2 and if ^n>3 are zero or small enough. The parameter region of applicability of our 
statement might be however even larger. The following reasons can be used in support 
of our statement. The task of finding time-periodic NLEs can be transformed into the 
task of solving a set of algebraic equations for the Fourier components A^i: Xi{t) = 
J2k Ai.iexp{iku}it) with k = 0, ±1, ±2, .... It was shown that in terms of certain separatrix 
manifolds the problem is generically solvable if the nonresonance condition is met [3]. 
Furthermore we have shown that the nonresonance condition can be obtained considering 



certain effective potentials of the original problem [4]. Finally we mention a recent NLE 
existence proof by MacKay and Aubry for the case of infinitesimal weak coupling ^{z), 
where an analytical continuation from the uncoupled case ^(z) =0 was carried out [5]. 

The time-periodic NLEs decay exponentially fast in space provided ^2 > [6]. They 
are exact solutions and thus never 'decay' in time without perturbations. As for per- 
turbations, we have shown that as long as the perturbations are not too strong, time- 
periodic NLEs can be locally perturbed into quasi-periodic NLEs (meaning that the 
time-dependence of a given variable is quasiperiodic). These quasiperiodic NLEs are 
not exact solutions of the equations of motion [6]. Quasiperiodic NLEs will radiate en- 
ergy and transform either into time-periodic NLEs or follow more complex scenarios [4]. 
Still quasiperiodic NLEs are meaningful objects, since their lifetimes can easily exceed 10^ 
periods of oscillations [4]. The effects of extended perturbations on time-periodic NLEs 
will be discussed in detail below. 

Time-periodic NLEs can be excited on any lattice site. Because of their robustness to 
perturbations these excitations can be easily observed in finite temperature simulations 
of example systems [1,7,8]. One could view NLEs as generalized quasiparticles with 
finite lifetime (at finite temperatures). On short time scales the quasiparticles behave as 
being uncoupled. On large time scales the coupling introduces relaxations, which lead to 
creation and annihilation of NLEs at different lattice sites. 

Recently Peyrard and coworkers have proposed a modified model [2,9], which is dif- 
ferent from (1),(2) in that the interaction potential $ contains terms {Xi -|-X;_i) in (2). 
Analysis of (2) yields then a remarkable change in the thermodynamic properties of H^ - 
the denaturation transition becomes exactly a first-order phase transition, whereas with- 
out these new terms the denaturation transition is only approximately of first order. As it 
emerges from the analysis [9], the change in the properties of the transition is sharp, even 
if the additional new terms are switched on perturbatively. Consequently the properties of 
single excitations will be only smoothly altered with variation of the perturbation. Thus 
at least for weak enough perturbations the NLE solutions will survive. Indeed numerical 
simulations [2] have shown that virtually no changes in the dynamical properties of the 
system appear as long as one does not consider the transition energy (or temperature) 
itself. Especially NLE structures have been observed both with or without the additional 
terms. 

Estimations for the parameters of our model can be found in the work of Prohofsky 
[10] and yield (l)2/v2 ^ 0.03. ..0.05. These parameter ranges certainly support the exis- 
tence of very localized NLEs [4]. In terms of the original physical interpretation of our 
model this means that DNA dynamics supports localized excitations of H-bonds between 
bases from a base pair. These excitations will have exponential decay in space (along 
the strands). By itself this statement does not imply any biological interpretation of our 
results. Of course we could speculate that an existing NLE might act as a precursor of 
a local DNA-opening. In fact it is even possible that NLEs exist where one H-bond is 
periodically opened and closed. Still the functioning of DNA is strongly controlled by 
certain enzymes. Consequently the most important question is whether it is possible to 
control the properties of the system in the controlled presence of NLEs. In other words, 
we think it might be much more reasonable to assume that enzymes interact with DNA 
strands and create NLEs for special purposes. If this hypothesis has some truth, then we 
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Fig. 1: Dependence of the onsite potential V on z (cf. Eq.(5)). 

indeed can assign biological functioning to an object which is described within nonlinear 
physics of many particles. 



3. FUNCTIONING WITH NLEs - PHONON SCATTERING 



We can roughly divide the possible functions an NLE fulfills into two classes - biological 
and physical. The first one is connected with specific chemical and biological processes 
and is certainly not accessible within our approach. What about the second possibility? 
Since our toy model is very simple, so are the possible properties of the system. One 
obvious property is energy or heat transport along the chains (strands). In the following 
we will show that the presence of an NLE can indeed control the heat flow and thus cre- 
ate a stationary nonequilibrium situation which is not accessible within a usual statistical 
approach (e.g. of finite temperature simulations). 

Let us confine ourselves without any limitation to a simple realization of model (2), 
such that the details of our results become more concrete. The onsite potential V{z) is 
choosen as 

' v^4z) = liiz - ly - ly , z<i 
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(5) 



To guarantee smoothness at 2; = 1 up to (including) third derivative we use < k < 0.25, 
a + K = 0.25, 2a P = 1. Variations of parameter k are not essential in the following, 
since we are not interested in the details of say the denaturation transition. The potential 
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Fig. 2: Numerical result of a time-periodic NLE solution with frequency coi = 1.3 for the model 
parameters as indicated in the text. The absolute values of the Fourier coefficients A^i of the 
periodic displacements Xi{t) are plotted on a logarithmic scale versus the lattice site I. The 
actual data are represented by open squares. Data for equal Fourier numbers k are connected 
with solid lines. The Fourier numbers from top to bottom are {k = 1, 0, 2, 3, 4, 5, ...). 



V{z) is shown in Fig.l for k = 0.2. The parameter V2 (cf. (3)) is equal to V2 = 2. 
We choose (f>2 = 0.1, 4>n^2 = in Eq.(3), so that we are in the parameter range as 
given in [10]. In that case numerical and analytical studies have revealed the existence of 
NLE solutions [4,11]. An example for a time-periodic NLE solution is shown in Fig. 2 by 
plotting the Fourier coefficients of every displacement variable Xi{t) versus lattice site. 
The logarithmic decay in space is clearly observed, in perfect agreement with analytical 
predictions [6]. These data were obtained using a specific discrete mapping which will be 
described elsewhere. Quasiperiodic NLE solutions have been also reported [4,11]. Their 
existence can be visualized with the help of Poincare maps [4]. The properties of NLEs 
can be studied with the help of reduced problems and effective potentials [4]. 

Let us turn to the question of heat flow control in the presence of a periodic NLE. We 
create a plane wave of phonons with a given wave number q and small amplitude. We 
want to measure the transmission coefficient of such a wave by a time-periodic NLE. Since 
the wave amplitude is assumed to be small, we can account for the scattering problem 
by adding small perturbations Si{t) to a time-periodic NLE solution Xi{t) and linearizing 
the equations of motion with respect to the perturbations: 
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Fig. 3: Discrete energy density distribution versus lattice site I of a scattering experiment after 
a waiting time of T = 12000. The infalling phonon wave has energy density e; = 10~^ and wave 
number q = 0.2n. The NLE is positioned at I = 1500. 

Let us first mention, that even if the nonresonance condition for the time-periodic NLE is 
met, v^e can not be sure that the NLE is stable in the presence of phonons. As a stability 
analysis v^ith respect to phonon perturbations shov^^s [4,11], stability of a time-periodic 
NLE is given if kuji ^ 2ujq. In the case of a stable NLE v^e expect the process of phonon 
scattering at an NLE to be elastic. The scattering of a phonon v^^ave (6) is then equivalent 
to the scattering of a single electron in a tight-binding model [12] v^ith a locahzed array of 
time-dependent diagonal defects. We do not knov^^ at present hov^ to analytically treat this 
problem of time-dependent multiple scattering. Consequently v^^e will present numerical 
data of the scattering. We use a time-periodic NLE v^^ith energy E = 0.256 and frequency 
Ui = 1.177. The phonon v^aves are created on one side of the NLE and have initial energy 
density 10^^ per lattice site. After an initial v^^aiting period (in order to make sure that 
phonon packets v^^ith other v^^ave number than the desired one travel away) we measure 
the transmitted wave over a time period of T = 12000. The system consists of 3000 lattice 
sites which is sufficiently large that we can exclude reflections from the boundary. The 
ratio of the transmitted to the incoming energy densities (squared transmission coefficient 
|tp) is then measured. In Fig. 3 we show the discrete energy density distribution ei as a 
function of the lattice site number / after waiting time T = 12000: 
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Here the NLE is positioned at lattice site I = 1500 and the incoming phonon wave 
(/ < 1500) with q = 0.27r travels to the right and is scattered by the NLE. First we find 
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Fig. 4: The squared absolute value of the transmittion coefficient |ip as a function of wave 
number q of the infalling phonon wave. The open triangles show the maxima (up) and minima 
(down) of |ip, the filled circles show the time-averaged value of |ip. The open squares correspond 
to the static linear Ersatzproblem. 

that the transmitted wave energy density is oscillating in time (as it should be because 
of the oscillating time-dependence of the diagonal defects in the scattering problem) (see 
Fig. 3, / > 1500). Secondly v^e find that the reflected v^ave does not loose its coherence 
v^ith the incoming wave - since we observe standing waves due to the interference between 
incoming and reflected waves (see Fig. 3, I < 1500). These standing wave structures are 
indeed stationary, as analogous pictures for different times clearly show. Note that the 
ordinate in Fig. 3 is plotted on a logarithmic scale. Since the corresponding wave length 
of the phonon is A = 10, one should expect the distance between two neighbouring knots 
of the standing wave (minima in the energy density) to be A/2 = 5. That is precisely 
observed in Fig. 3 (/ < 1500). 

But the most important information is the dependence of the squared transmission 
coefficient on the wave number of the incoming wave as shown in Fig. 4. We observe a 
dramatic decrease of |tp over orders of magnitude with increasing wave number q. Only 
for very small wave numbers (large wavelength) do we find appreciable amounts of trans- 
mitted wave intensity. This result might appear to be obvious - since large wavelength 
phonons do not 'feel' the finite size perturbation by the NLE. This is however wrong, as 
the simple example of one static diagonal defect shows - the transmittion coefficient for 
q = is zero as well as for g = tt (because of the singularities in the phonon density of 
states) [12]. 

If we replace the time-dependent defects in equation (6) by their time-averaged val- 
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ues, we obtain a linear Ersatzproblem. The numerical investigation of the phonon scat- 
tering for this Ersatzproblem yields the open squares in Fig. 4. Allthough the qualita- 
tive g-dependence seems to be reproduced, the quantitative difference between the time- 
dependent and the time-independent scattering problems is huge. Consequently we are 
confronted with a rather subtle result of time-dependent multiple scattering events, which 
leads, as the numerical result clearly shows, to the fact that phonons with almost any wave 
number are strongly reflected by a time-periodic NLE. Thus we can prevent the heat flux 
from penetrating certain areas of our chain by creating two NLEs on the borders of de- 
sired region or conversely we can trap energy in a region between two NLEs.. This NLE 
creation does not require an overall rearrangement of the chain as in the case of a kink 
but only a local excitation of one or a few degrees of freedom is needed. 

If we consider larger interaction ranges than nearest neighbour interaction, we do not 
expect that the NLE properties will change drastically, as long as we have finite range 
interaction. As analytical considerations indicate [3,6], there is no basic difference between 
say next nearest neighbour interaction and nearest neighbour interaction with respect to 
the existence of NLEs. 

Let us return to the problem of DNA functioning. One possible role for NLEs in DNA 
functioning could be an enzyme which could create NLEs in order to prevent certain ar- 
eas of DNA from being penetrated by heat flux or to store energy in certain areas. This 
could be an example of functioning of DNA describable on the level of nonlinear physics 
of many particles. Of course there can be many more possibilities for functioning using 
NLEs, however to describe these we would have to leave the abstract level of nonlinear 
physics and become much more specific in our understanding of the system. It makes no 
sense to speculate about these possibilities without experimental evidence. 
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